arXiv: 1509.04003v2 [quant-ph] 4 Aug 2016 


Ultra-small time-delay estimation via weak measurement technique with postselection 


Chen Fang*/ Jing-Zheng Huang*/' Yang Yu,' Qinzheng Li,' and Guihua Zeng^'’^ 

^ State Key Laboratory of Advanced Optical Communication Systems and Networks, 

Shanghai Key Laboratory on Navigation and Location-based Service, 
and Center of Quantum Information Sensing and Processing, 

Shanghai Jiao Tong University, Shanghai 200240, China 
^College of Information Science and Technology, Northwest University, Xian 7I0I27, Shaanxi, China 

(Dated: August 5, 2016) 

Weak measurement is a novel technique for parameter estimation with higher precision. In this paper we 
develop a general theory for the parameter estimation based on weak measurement technique with arbitrary 
postselection. The previous weak value amplification model and the joint weak measurement model are two 
special cases in our theory. Applying the developed theory, the time-delay estimation is investigated in both 
theory and experiment. Experimental results shows that when the time-delay is ultra small, the joint weak 
measurement scheme outperforms the weak value amplification scheme, and is robust against not only the 
misalignment errors but also the wavelength-dependence of the optical components. These results are consistent 
with the theoretical predictions that has not been verified by any experiment before. 

PACS numbers: 


I. INTRODUCTION 

As an advanced technique that provides higher sensitivity 
for parameter estimation, the weak measurement scheme has 
attracted much attention in recent years iQH . Commonly, the 
weak measurement technique with postselection involves two 
physical systems, called ancillary system and pointer. These 
systems are weakly interacted which may be described by a 
coupling parameter. To estimate this extremely small parame¬ 
ter with higher sensitivity, the ancillary system is postselected 
by two orthogonal states, and a followed measurement oper¬ 
ated on the pointer system provides available information for 
estimating the coupling parameter. 

By far, two special weak measurement techniques with 
postselection, i.e., the weak-value amplification (WVA) iQl] 
and the joint weak measurement (JWM) H, have been in¬ 
vestigated. The main characteristic of these techniques is that 
the postselection can improve the estimation precision in the 
presence of technical restrictions In the WVA scheme, 

the measurement result is recorded when the ancillary sys¬ 
tem is successfully postselected by one of the postselected 
states, which is chosen almost orthogonal to the initial state. 
This way leads the postselection probability to be extremely 
low. Using the recorded results, the coupling parameter is es¬ 
timated by simply averaging the measurement results ifl^ . In¬ 
vestigation shows that the WVA is useful in suppressing some 
technical noises and realistic limitations, such as detector sat¬ 
uration, correlated noises and angular beam jitter i^. Subse¬ 
quently, it is easier to reach the quantum standard limit with 
common experimental equipments in the WVA scheme ifnl] . 
Moreover, in the scenarios of phase-shift and time-delay mea¬ 
surements, the WVA technique may provide at least 3-orders 
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of magnitude improvement on estimation precision over the 
standard interferometr y lUOl] when similar alignment errors 
are taken into account lll^fl^ . Of course, there are also two 
disadvantages, i.e., the restrictions on the estimati on p recision 
due to systematic errors and alignment errors ifiolll^ . and the 
limitations on the efficiency due to the waste of large amount 
of unselected events after postselection operations lU^ . 

To overcome the disadvantages in the WVA scheme, a nat¬ 
ural way is to collect the unselected data iflsilT^ . This leads 
the so-called JWM scheme |01, in which all measurement re¬ 
sults are recorded when the events nearly equiprobably post- 
selected by the two orthogonal postselected states. With the 
measurement results, the coupling parameter is estimated em¬ 
ploying the maximum likelihood estimation method. Ref.|01 
proved that by carrying out the JWM on all events, one can not 
only increase the efficiency but also remove systematic errors 
and alignment errors. E^ecially, the JWM technique has no 
ultimate precision limit lalSDj which is very useful in practice. 

In this paper, we develop a general theory of weak mea¬ 
surement technique with arbitrary postselection. Since the in¬ 
volved postselections in the WVA scheme and in the JWM 
scheme are extremely unbalanced and almost balanced, re¬ 
spectively, they are actually two special cases in our general 
model.Our work is distinguished from the previous works that 
concerned in the WVA regime (13, [Hi- With the proposed 
general theory, the ultra-small time-delay estimation is inves¬ 
tigated. Especially, a novel experiment for time-delay esti¬ 
mation making use of the weak measurement technique with 
postselection is firstly presented 113] . Our experiment verifies 
exactly the theoretical predictions in Ref.i^, and we find that 
when the time-delay is ultra small, the JWM technique out¬ 
performs the WVA technique, and is robust against not only 
the misalignment errors but also the wavelength-dependence 
of the optical components. 

The paper is structured as follows. In Secjlll a general 
theory for parameter estimation based on weak measurement 
technique with arbitrary postselection is proposed and then is 
applied in the time-delay estimation. In Secln!] an ultra-small 
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time-delay estimation experiment is presented. Finally, the 
conclusions are drawn in Sec IIVI 


II. THEORY 


In this section, we develop a general theory for parameter 
estimation based on weak measurement technique with arbi¬ 
trary postselection. Then we apply it in the ultra-small time- 
delay estimation. 


A. General Theory of Weak-Measurement-Based Parameter 
Estimation 

Consider a weak measurement scenario involving a two- 
level ancillary system prepared in state \ifi) and a pointer with 
a continuous degree of freedom in state |0) = f dp4>(p)jp}, 
where p is a continuous variable and (/)(p) is the corresponding 
wave function. The interaction between the system and the 
pointer is described by an unitary operator 

U.nt = (1) 

where A acts on the system with eigenvalues of 1 and — l,p 
acts on the pointer, and g is the coupling parameter which in¬ 
dicates the coupling strength. This operator can be expressed 
as JUl]: 


Uint = e = cos{gp)i - i sm{gp)A. (2) 

After the interaction, an initial product state |'I') = |<pi) 0 
I (j)) evolves to 

I'k') = Uint\Pi)\(l>) ^ ^2) 

= [cos{gp)\(pi) -ism{gp)A\(pi)]\(j)). 

The ancillary system is then postselected by two orthogonal 
states, i.e., |p/i) and lp/ 2 ), which converts the pointer state 
to 


= (<F/jl<Fz)[cos(pp) - isin{gp)A^j]\(j)), 


(4) 


where Qj{p){j = 1, 2) is the probability distribution of ob¬ 
serving p in the experiment, which in principle converges to 
Pj (p) as the number of measured events is increased. Then, 
the maximum likelihood estimate of g can be achieved by 
solving the likelihood equation dL{g)/dg = 0 ll2^ . In our 
case, it yields 


— V' r p[sin2(gp)(|A,^jF-l)-r2cos2(gp)ImA„3-] 

J PWj \P ) cos^ {gp)-\-sin^ {gp) \ Aw j\^ -\-sin 2{gp)lm Aw j 

= 0 . 

(7) 

As the weak measurement technique requires pp <C 1 and 
\Aw\gp ^ 1 ll^ . Eq.(|7]i can be simplified as 


Ej / dpQjip) 


-2lmA^jP^ g'^ + {\A^j\^-l)p'^ g+lmA^jP 
(\An,j\^-l)p^g^+21mAu,jgp+l 


= 0 . 


( 8 ) 


Interestingly, Eq.® shows that prior knowledge about Po{p) 
is unnecessary. 

Generally, the analytical solutions of Eq.® is difficult to 
obtain. Eortunately, under the condition of \Auj\gp <C 1, 
Eq.® can be approximately simplified to a quartic equation 


Ag^ + Bg^ + Cg^ + Dg + E = 0, (9) 


where 

A = J2j I dpQj{p)p^[2{\Ayjj\'^ - l)ImA™j], 

^ = T,g IdpQj{p)p'^AA^A^j)'^ - - 1 ) 2 ], 

C = J2j I dpQi(p)p^[ImAu,j(l - 3\A^j\% 

dd = J2j I dpQj{p)p'^[{\A^j\'^ - 1) - 2ImA^,j], and 
E = J2j f dpQj{p)plinA^j. 

As g is extremely small, one may obtain the first-order ap¬ 
proximate solution 

„(i) = _E =__ .,01 

D E[2(W^,)2-|A^,P + l](p2)/ 

where (p)j = f dpQj{p)p and (p^)^ = f dpQj{p)p^. 
Eg. ([Tol l provides an appropriate estimate of g with very high 
precision. 


B. Time-delay estimation 


where 7 = 1,2, and Ayji = 

weak value yj]. Consequently, the probability distribution of 
p associated with the pointer becomes 

pap) ^\m)? ... 

= \lpf3\Pi)?Po{p)Q{P^9), 

where PAp) = the initial probability distribution 

of the pointer, and Cj{p,g) = cos^(gp) + siv?{gp)\Awj\^ -f 
sin 2 (( 7 p)ImAujj is the part changes the shape of the initial 
probability distribution. 

To achieve an unbias estimate of g, we apply the maximum- 
likelihood estimation method llEI] . Eirst, we construct the log- 
likelihood estimator as 

Pg) =T,j I dpQjip) log Pjip), (6) 


As an application of our developed theory, we consider the 
time-delay estimation described in Eig[T] In this scheme, an 
ultra small time delay of r between two orthogonal polariza¬ 
tions, which are denoted as horizontal (H) and vertical (V), 
is involved by a thin birefringent crystal. The corresponding 
unitary operator is given by 

( 11 ) 

where az = \H){H\ — |E)(V| is the Z operator acting on 
the polarization states and uj represents the frequencies of the 
photons. In this scenario the parameters g and p in the general 
model are replaced by time delay r and frequency w, respec¬ 
tively, and the photon polarizations and frequencies are cho¬ 
sen as the ancillary system variable and the pointer variable, 
respectively. 
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optical axis 



FIG. 1; (Color online). Time-delay estimation. Photons emitted from 
a light source is modulated to be a linear superposition state of two 
orthogonal polarizations, H and V, by a linear polarizer. After pass¬ 
ing through a thin birefringent crystal, a small time delay is involved 
between the two polarizations. The followed circular polarization 
beam splitter(PBS) postselects two polarization states, and two spec¬ 
trometers (D1 and D2) are placed to measure the spectrum of the 
output lights. 





FIG. 2: (Color online). Representation on the Bloch sphere of the 
pre-selected and post-selected states, where the subscripts WVA and 
JWM stand for two schemes. 


Given an initial state of and post-selected 

states of \ipfi) = and \ipf 2 ) = ’ we get 

the weak values 


A _ _ l-e-'* 

A _ (yf2|A|yi) _ 1+e-^'t’ 


i tan(^) 
—i cot(|). 


( 12 ) 


The postselection probabilities are given by 

Pfi =\{^fiWi)? S du:Po{uj)C,i{u:,T) 

= cos^ 2 I — + 2tan|a;oT]; 

Pf2 =\{<Pf2\‘Pi)\'^ J dujPo{uj)C2i^,T) 

= sin^ |[1 + (cot^ I - l)woT^ - 2 cot |wor], 

(13) 

where ujq = f Po{u!)ujduj represents the initial average fre¬ 
quency of light before the interaction. In practice, Pji and 
Pf 2 can be estimated from experimental data by f Qi(aj)duj 
and f Q 2 (cu)duj, respectively. 

Inserting Ea. (fT^ into Eg. lfTOl l gives 

(1) _ sin^ j(i^)i-cos^ |(a;)2 ,, ,, 

tan ^ ^ ' 

We note that estimating r with Ea. (fl4l l requires prior knowl¬ 
edge on (j), which is, however, not necessary in high preci¬ 
sion for achieving high signal-to-noise ratiol^. Eurther dis¬ 
cussions can be found in Sec.II.C. 

Now we consider applications of the proposed theory 
model in two specific schemes, i.e., the JWM scheme and the 
WVA scheme. In the JWM scheme, the light intensities de¬ 
tected by two spectrometers are nearly equal, which requires 
(p Ri f, the pre-selected and post-selected states can be seen 
in EiglU The postselection probabilities of two output ports 
are 


In the JWM scheme, the weak value is not anomalous, thus 
one can not observe an significant shift on the pointer’s wave 
function, as it is the case in WVA. In fact, the measured shift 
is existed in the probability difference ll^ . 

Then, Ea. (fT4l i can be simplified as 

_ ^ P.fi(^)2-f.f2<‘^)i ('If;') 

Ijw ~ 

Employing Ea. ifThl l the time-delay r can be estimated without 
being precisely aware of p. We note that this result is conflict 
with the formula given in i^, which is 


^ I Auj)2 - (cu)l _ Pf2 - Pfl , 

4^ Aa;2 Aw 

In SeclIII]we will further verify the correctness of our formula 
by experiment. 

Next we consider the application of the proposed theory 
model in WVA technique. In the WVA scheme, most of the 
light output from port 1, which requires ^ r; 0 as shown in 
EiglU The postselection probabilities of the two output ports 
are 


Pfl ~ (1 - 

Pf2 r:^[1-^Wot]. 

In this case, Ea. (fT4l) gives 


(17) 


'^wva 


{cj^)^ + {ui^)2-2uio{uj)2 ’ 


( 18 ) 


where 


and 


r- f Qj{uj)ujduj 


{'P)j/Pfh 


Pfl « cos^ |(1-F 2tan^woT). 
P /2 ~ sin^ |(1 — 2 cot 


-r^ ^ J Qjiuj)uj^duj 

^ f Qj{^)duj 


(15) 


)j/Pfj- 
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In Eg.lfTsTl. we have (w)^ — {uj)^ ~ wq — 

+ (a ;^)2 — 2wq{ui)2 ~ 2(Aa;^). Reminding that = 
-tcot(|) Ri we get 

Suj Ri 2Test/mAu,2(Aw^), (19) 

which is the familiar formula derived in the weak value am¬ 
plification scheme ifl^ . 

III. EXPERIMENT 

In this section, we perform an experiment to demonstrate 
the time delay estimation based on the theory model presented 
in Sec III Bl Making use of the experimental results, we study 
the performances of the JWM and WVA schemes, so that we 
can verify the theoretical predictions about the advantages of 
JWM presented in Ref.ll^. 

A. Experimental Setup 



FIG. 3: (Color online). Experimental setup of the time-delay es¬ 
timation. Photons emitted from a light emitted diode (LED) with 
center wavelength of 780nm and line width of AA = 17.6nm are 
collimated hy an aspheric lens and a convex plane lens then en¬ 
ter the Glan-Taylor prism, which preselects the polarization state 
\ipi) — -^{\H) + \V)). After that, the light passes through two com¬ 
pound binary zero-order half-wave plates(HWPs) with their optical 
axes(OAs) perpendicular to each other and at 45° to the polarization 
of the prism.The second HWP is pivoted by a small angle 9 around its 
OA. Postselction is realized by a combination of a compound binary 
zero-order quarter-wave plate(QWP) and a polarization beam split- 
ter(PBS).The photons are split into two orthogonal polarized beams 
and focused into optical fibers by lens, then sent to two spectrometers 
with sampling resolution of 0.1 nm and range of 690-900 nm. 

The experimental setup is described in Figl3 in which a 
thin birefringent crystal is used to induce the time delay. We 
note that it is very difficult to produce and manipulate this 
kind of plate in practice. Follow the idea in Ref. lU^ . we use a 
double-plate system with placing their optical axes(OAs) per¬ 
pendicularly to each other to equivalently realize this effect. 
Restricted by the laboratory condition, we use the binary com¬ 
pound zero-order half-wave plates (HWPs) in our experiment. 


These plates are built by combing two multi-order wave plates 
and aligning the fast axis of one plate with the slow axis of the 
other to obtain the zero-order phase delay. In contrast to the 
true zero-order plates used in 111311 . they are insensitive to the 
temperature changes. We place both plates perpendicularly to 
the light velocity to cancel their phase delay. The delay is re¬ 
alized by pivoting the second plate around its OA by a tiny 
angle 9. This pivot increases the optical path of the second 
HWP, which makes the double-HWPs system equivalent to a 
very thin plate of the same material with the OA orienting as 
the tilted one. The relationship between the time delay t and 
the tilt angle 6 is 

-(ne-no)/te^ 

where Ue, rio and n are the refractive indexes of quartz for 
extraordinary light, ordinary light and average, respectively, h 
is the thickness of the plate, c is the light speed in vacuum and 
A is the wavelength of the light. For detailed derivations, see 
Appendix. A and Ref, . 

Finally, in order to finish the postselection, the photons se¬ 
quentially enter a quarter-wave plate (QWP) and a polariza¬ 
tion beam splitter (PBS). Originally, the OA of the QWP is 
perpendicular to the axis of the Glan-Tylor prism and the PBS 
splits the photons to polarization states \H) and \V). And 
then, the PBS is rotated by an angle of 7 = 7 r /4 — ^ ac¬ 
cording to the GT prism’s OA orientation to postselect fi¬ 
nal states of |(p/i) = C 0 S 7 |iT) -F sin 7 | 17 ) and |<p/ 2 ) = 
sin7|iT) — cos7|y). 

In our setting, when (f) m the light intensities go out from 
the two output ports are nearly the same, which corresponds 
to the JWM proposal i^. In contrast when (jj m 0, almost all 
of the light go out from one output port, which corresponds 
to the weak-value amplification scheme. In the following, we 
will present the experimental results of these two schemes to 
verify the theory derived in SecHIl 

B. Joint Weak Measurement v.s. Weak Value Amplification 

Firstly, we consider the JWM technique where the post¬ 
selection is almost balanced. In this scheme, we set (p Ri 
-F 0.071) rad (by rotating the PBS with an angle of 
—0.071 rad), and the postselection probabilities for each out¬ 
put port can be calculated by Fq.lfTSll. In Fig|4] according to 
the tilted angle we made, different estimations about the time- 
delay are shown, where the theoretical predictions drawn by 
blue line are calculated by Fa.(l20ll. Among these estimations, 
blue triangles are the complete solutions to Fq.®, when red 
circles and black crosses are the approximate solutions given 
by Fa. (fl4l i and Fa.ifThll respectively. Compared to Fa. (fT4l i. 
Fg. lfThl l does not require the prior knowledge about p, this pa¬ 
rameter is estimated by the light intensity difference of the 
two output ports. However Fg. lfThb provides an estimation 
with less accuracy, because the environment effects, such as 
the fluctuation of the temperature, make the output intensity 
of the LFD unstable . 

Then, we consider the WVA technique where the postselec¬ 
tion is extremely unbalanced. In this scheme, we set (p Ri 0.03 
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FIG. 4: (Color online.) Relation of 6 and r in the joint weak mea¬ 
surement scheme. Blue line presents the theoretical prediction from 
Eg. OOt . blue triangles, red circles and black crosses are estimates of 
T from experimental data, calculated by the most rigorous equation 
[Eq.®], the first-order approximation lEa.l ll4H and the simplified 
approximation lEg.JlfiH respectively. 


(by rotating the PBS for 7r/4 — 0.03 rad), and the postselec¬ 
tion probabilities for each output port, which can be calculated 
by Ea.(fT7ll. are extremely unbalanced. In Fig® according to 
the tilted angle we made different estimations about the time- 
delay are shown, where the theoretical predictions drawn by 
blue line are calculated by Ea. (l20l i. Among these estimations, 
blue triangles are the complete solutions to Eq.®, when red 
circles are the approximate solutions given by Ea. (fT4l) . 

We can find that the red circles fit well with the theory when 
T is less than 0.5 x 10“^^ but the deviation increases along 
with the growing r, because the effect of wavelength depen¬ 
dency of the quarter-wave plate was not taken into account in 
Eg. lfTSl) (see Appendix. B for details). For more discussions 
about this phenomenon, we refer the readers to Refs.J^^-H^. 


C. Discussions 

Utilizing the imaginary part of weak value provides a large 
amplification factor for parameter estimation. For this pur¬ 
pose, a quarter-wave plate for adding a circular component 
on the polarization is required. However, as is pointed out 
in 0, the quarter-wave plate involves large uncertainty to 
the polarization that is difficult to compensate, because of the 
wide spectrum of the light source. 

There are two possible ways to solve this problem: get rid 
of the quarter-wave plate with the price of much smaller weak 
value, or apply joint weak measurement method to suppress 
the wavelength-dependence effect. 


FIG. 5: (Color online.) Relation of 9 and r in the weak value am¬ 
plification scheme. Blue line presents the theoretical prediction from 
Eg. 1201 . blue triangles and red circles are the estimates of r from 
experimental data, calculated by the most rigorous equation [Eq.®] 
and the first-order approximation lEg.l lldH respectively. 

In our experimental setting (see Fig®, removing the 
quarter-wave plate reverts a similar setup as that in ifT^ . Al¬ 
ternatively, we can collect the data from both output sides and 
apply Ea.tfTOli for parameter estimation. When the condition 
\A^\ujot <C 1 is satisfied, the estimation fits quite well to the 
theoretical prediction (see Fig®. 

In ifT^ . the uncertainty of measuring a (equals uiifr in our 
paper) depends on the spectrometer resolution A((5A) and the 
postselection parameter ^ (equals 27 in our paper). The un¬ 
certainty of a, which is denoted as Aa, can be derived as lfl^ : 

A«=4^-^^^-A(a). (21) 

According to this formula, the optimal precision can be 
achieved when /3 ~ a, but it is very impractical because the 
value of a is unknown before it is measured. Alternatively, 
the authors pointed out that for measuring small ol one can set 
/3 = 0 and obtain Aa ~ 0.1a when A((5A) = O.lnm. How¬ 
ever, an important factor, the misalignment of /3, had not been 
considered. In fact, it sets a limitation on the minimum de¬ 
tectable phase when using weak value amplification lfl^ . As¬ 
sume that the actual value of /3 is e when we set it to be 0, and 
consider amm = Aa as the minimum detectable value of the 
measured phase, we can derived from Ea. (l2Tl l that 

_ / vu:4C-2C+r (22) 

where we define C = A((5A) for simplicity. By setting 

A{SX) = O.lnm, which is assumed in lUAII and practically 
applied in our experiment, we get amin — 3.7e. 

On the other hand, by using the joint weak measurement 
scheme proposed in |@] and demonstrated in Sec lIIIBI the un- 
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certainty of measured phase is insensitive to the alignment er¬ 
ror. In our experiment, we set (jj = 0.03. The exact value of 
(p can be estimated by maximizing Eq.®, but it is not neces¬ 
sary to do so for achieving high precision. In Figl6] we cal¬ 
culate the signal-to-noise ratio(SNR) when estimating r by 
Eg. lfTOl) and setting (p = 0.03, 0.05 and 0.08 rad respectively. 
The SNR drops slightly even when our prediction about (p has 
a deviation of 0.05 rad, and remains higher than 10 dB for 
a > 0.002. It shows that the estimation results are insensitive 
to (p in balanced postselection scheme. In contrast, for weak 
value amplification in unbalanced postselection scheme, for a 
typical deviation of e = 0.0027 11311 , the SNR drops down to 0 
dB when a = 0.01, much worse than applying the joint weak 
measurement method. 
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FIG. 6: Signal-to-noise ratio in the time-delay estimation experiment 
with joint weak measurement. Red circles; p = 0.03; Blue stars: 
p = 0.05; Black crosses: p = 0.08. 

As is predicted by theory|@] and presented in Fig|6l the 
SNR of joint weak measurement scheme remains stable when 
a grows, because the estimation error increases along with the 
signal. As in the weak value amplification scheme the SNR in¬ 
creases along with the signal 11^ . it implies that at some point 
the weak value amplification scheme will achieve higher SNR 
than the joint weak measurement scheme. In that case, one 
can adopt the weak value amplification scheme for better per¬ 
formance, however the harm of the QWP must be overcome 
by, e.g., applying a broadband wave plate. This is beyond the 
scope of our current work. 

Finally, it is useful to compared the JWM and WVA with 
other schemes for practical purpose. In Ref. m the authors 
compared the WVA and standard interferometric scheme for 
phase shift estimation, and proved that when similar align¬ 
ment errors are taken into account, the estimation precision of 
WVA is at least 3-orders of magnitude higher than the stan¬ 
dard interferometric scheme. Combining with the above dis¬ 
cussions, it is reasonable to expect that JWM has the highest 
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precision for estimating extremely small time delay among 
these three schemes. 


IV. CONCLUSION 

In summary, we present a general theory for parameter es¬ 
timation based on weak measurement with arbitrary postse¬ 
lection. Applying this theory, we study the time-delay esti¬ 
mation in both theory and experiment, with especial interest 
in two specific schemes, i.e., the WVA scheme and the JWM 
scheme. Through the experimental results, we find that the 
JWM scheme outperforms the WVA scheme when the time- 
delay is ultra small. These results support the theoretical pre¬ 
dictions presented in Ref.|@] that has not been verified by 
any experiment before. Furthermore, the JWM technique is 
robust against not only the misalignment errors but also the 
wavelength-dependence of the optical components, which ap¬ 
pears to be an advanced parameter estimation approach for 
achieving higher precision with convenient laboratory equip¬ 
ments. 
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Appendix A: Phase retardation of compound binary zero-order 
wave-plate 

When light propagates through a uniaxial crystal plate, the 
phase delay between the ordinary light and extraordinary light 
is. 


S = ^{ne — no)L sin^ p, 

where A is the wavelength,no and Ue stand for refractive index 
of ordinary and extraordinary light,respectively.L = 
means the path length propagating in the birefringent crys¬ 
tal. f] is the angle between the normal line and the OA,iji 
is the angle between the light velocity in the crystal and the 
normal line.And h is the thickness of the plate. For sim¬ 
plicity, we establish a Cartesian coordinates showed in FigjT] 
to calculate C = L sin^ rj in the crystal which we call the 
effective path length in the next. The origin is coincident 
with the point of incidence.X-Z plane is on the surface of the 
plate.And the OA is paralleled with Z-axis. It’s obvious that 
sin^ t] = + y^,and 


C = L sin^ p 


2 I 2 

X +y 


(A2) 


The exit point is(a;, h, z).And the effective path length could 
be obtain by replacing h by y.In our experiment, the pivot is 
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FIG. 7: (Color online). Coordinate system of single birefringent 
waveplate. 


around the OA, we define the azimuth angle and elevation 
angle as ^ and ip respectively, so 


^2 _|_ /j2 _ j (A3a) 

^2 _|_ /j2 _ 1^2 j (,Qg2 


and we could deduce that 


x'^ +y‘^ + z 


2 _ sin^ ^sin^'!/>) 

cos^~^~cos^~2pJ) ’ 


(A3b) 


inserting Ea lA2l and Ea lA3l into Ea lAll we could get the gen¬ 
eral expression of the phase retardance of oblique incidence 
on a crystal plate 


6 =^{ne- no)C 


2'K{ne — rLo)h cos ■0 
cos ^(1 —sin^ ^ sin^ -i/i) 


(A4) 


The compound binary plate is built by combining two 
multi-order wave plates to obtain an zero-order phase delay 
by aligning the fast axis of one plate with the slow axis of the 
other.The retardance could be expressed as 

6 = ^{ne-no){hi-h2). (A5) 


90°. Let the subscripts denote the two plates. We neglect the 
refringence between the two plates,and suppose that ^ in EiglS] 
is approximately equal to ip in Eig|7l we can derive that 


Zi/hi = X2/h2, 
xi/hi = -22/^2- 


(A6) 


The phase delay of the compound plate is 

S=^{ne-no){Cl-C2), (A7) 


where 


Cl = 
C2 = 


xl+hf 

yCf+iif+zf ’ 
h2{xl+hl) 

hly^xl+hl+zf ’ 


(A8) 


combining with Ea lA2l and taking account of the relation¬ 
ship of those angles, we found the general expression of the 
phase retardance of the compound plate; 


j _ 27r(ne —no)(/i-i COS t/j/COS ^ —/i2 COS COS-i/j) (A9) 

when Ip = 0 and ^ is small, Ea lA9l becomes 


<5 


27r(ne —no) 

27r(ne —no) 
A 


(/ll/ COS^ — /12 COS^) 

\{hi — ( 12 ) + + ^ 2 ) 72 ], 


(AlOa) 


and when ^ = 0 and ip is small, Ea lA9l at)r)roximates to 


S 


27r(n^ cosip — h2/cOS1p) 

Pi _ h,) - + h2)/2]. 


(AlOb) 


In Eq IAlOl the 1st term is the plate’s retardance of vertical 
incidence,the 2nd term means the phase delay of oblique inci¬ 
dence.Considering the relationship between the incidence an¬ 
gle 9, Ip and we could rewrite the 2nd term as 



FIG. 8: (Color online). Coordinate system of the second plate in a 
compound wave plate. 

FigUl represent the Cartesian coordinates when light go 
through the second plate while the EiglT] still shows the inci¬ 
dence on the first plate.The origin overlaps with the exit point 
of the first plate.Since the OA of the second plate is perpen¬ 
dicular to the first one’s, the coordinate system is rotating by 


_ ±Tr{ne-no)(hi+h2)6^ (All) 

when '0 = 0 it has the positive value and the negative value 
could be obtain while ^ = 0. In the range of our LED spec¬ 
trum ,ne — Uo is almost constant, in other words, the optical 
length in the crystal is unsensitive to the variation of wave¬ 
length. With the cancelling effect of the two waveplates in our 
experiment, the time delay induced by the pivot is 

_A 5 A _ ±{ne-na)(hi+h2)0^ A 1 "rt 

-271115- 

where c is the the speed of light in vacuum. 

Appendix B: Effect of wavelength-dependent quarter-wave plate 

The postselection in our experiment consists of a quarter- 
wave plate and a linear polarization beam splitter(PBS). Eor 
using ideal quarter-wave plate Ea. dTSl i is established, however 
in practice optical components are wavelength-dependent. A 
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quarter-wave plate with its optical axis at 45' direction and 
central angular frequency of ujq can be described by a Jones 
matrix as; 


t^A/4 2 yi I J y 0 

_ / coswro —isinwTo \ 
I —isinwTo coswtq J 


1 

-1 



(Bl) 


where tq — jjuJo, and w is the angular frequency of the light. 

Denote the linear polarization state selected by PBS is \ipf), 
the postselection of |(^^) can be divided as 


= {{^f\U^/4)U^nt\^^M■ (B2) 


and 


\^'f2) = (cos wtq sin7 — i sinwTo cos7)|i?) 

-I- (—coswtq C0S7-f i sinwTo sin7)|14). 


Consequently, the weak values can be calculated as 


_ (cos cDTo+i sin ojTo)(cos 7—sin 7) 

(cos LOTQ — i sin a;ro)(cos 7+sin 7) 
_ COS 7 —sin7 22ajrn 
cos 7+sin 7 


(B3) 


Obviously, |+) = [/|^^|(/?/) is frequency-dependent. Com¬ 
bining with Eg. dBll i and the following expressions, i.e., 

1‘P/i) = cos7|JJ) -I- sin7|C) 


A _ 

^“2 - 

_ (cosojTo+z sincjTo)(cos 7+sin7) 

(cos lotq — i sin a;ro)(sin 7—cos 7) 
_ cos7+sin7 ^22ajrn 
sin 7—cos 7 


(B4) 


and 

\^Pf2) = sin 7 | 7 J) - COS7+), 

we can derive 

1+1) = (coswTo COS7-I-isinwTo sin7)|iJ) 
-f (cos wtq sin 7 -f i sin ujtq cos 7) +). 


In this case, A^j can not be taken outside of the integral of 
Eq.®. Compare the results shown in Eig|4]and EiglSjwe can 
see that in the balanced postselection scheme, Ea. ifTOl) pro¬ 
vides a good approximation, while in the unbalanced post¬ 
selection scheme the wavelength-dependency of the quarter- 
wave plate involves an unnegligible deviation. 
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